UMN-TH-2302/04 
TPI-MINN-04/13 
astro-ph/0404257 
April 2004 



The Fine-structure Constant as a Probe of Chemical Evolution 
and AGB Nucleosynthesis in Damped Lyman-a Systems 

Timothy P. Ashenfelter, Grant J. Mathews 

Department of Physics, Center for Astrophysics (CANDU) and Joint Institute for Nuclear 

Astrophysics (JIN A), 
- University of Notre Dame, Notre Dame, IN, 46556 

o : 

O ■ Keith A. Ohve 

(N ; 

! William I. Fine Theoretical Physics Institute, 

1—5 ! University of Minnesota, Minneapolis, MN 55455, USA 

i> : 

(SI : ABSTRACT 
> ■ 
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lO . Evidence from a large sample of quasar absorption-line spectra in damped 

^ ! Lyman-a systems has suggested a possible time variation of the fine structure 

O ! constant a. The most statistically significant portion of this sample involves 

Q ! the comparison of Mg and Fe wavelength shifts using the many-multiplet (MM) 

^ ! method. However, the sensitivity of this method to the abundance of heavy 

^! isotopes, especially Mg, is enough to imitate an apparent variation in a in the 

2 ! redshift range 0.5 < z < 1.8. We implement recent yields of intermediate mass 

^ ! (IM) stars into a chemical evolution model and show that the ensuing isotope 

! distribution of Mg can account for the observed variation in a provided the early 

• ^ ! IMF was particularly rich in intermediate mass stars (or the heavy Mg isotope 

^ ! yields from AGB stars are even higher than in present-day models). As such, 

i these observations of quasar absorption spectra can be used to probe the nucle- 

osynthetic history of low-metallicity damped Lyman-a systems in the redshift 
range 0.5 < z < 1.8. This analysis, in conjunction with other abundance mea- 
surements of low-metallicity systems, reinforces the mounting evidence that star 
formation at low metallicities may have been strongly influenced by a population 
of IM stars. Such IM stars have a significant influence on other abundances, par- 
ticularly nitrogen. We constrain our models with independent measurements of 
N, Si, and Fe in damped Lyman-alpha systems as well as C/0 in low-metallicity 
stars. In this way, we obtain consistent model parameters for this chemical- 
evolution interpretation of the MM method results. 



Subject headings: cosmological parameters — - galaxies: abundances — galaxies: 
evolution — stars: AGB — quasars: absorption lines 
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1. Introduction 

The origin and dynamics of the fundamental constants of Nature is one of the deepest 
questions in physics. One of the most widely held tenets in physics is that the laws of nature 
are universal, constant, and favor symmetries. Nevertheless, in various unified theories 
(including string theory), gauge and Yukawa couplings often appear as dynamical variables 
which are only "fixed" when a related scalar field (such as a dilaton in string theory) picks up 
a vacuum expectation value. While one may naturally expect that couplings become constant 
at or near the unification scale, only experimental evidence can constrain the degree to which 
these constants vary at late times. 

In this context there has been considerable excitement in recent years over the prospect 
that a time variation in the dimensionless fine structure constant, a, may have been observed 
(Webb et al. 1999, Murphy et al. 2001a,b, Murphy et al. 2003a). This evidence is based upon 
an application of the "many multiplet" (MM) method to quasar absorption lines in damped 
Lyman-a systems (DLAs). Attempts at constraining or measuring time variations of a in 
quasar absorption-line spectra have a long history going back to work by Bahcall & Salpeter 
(1965) using 01 1 1 and Nell I emission lines. This method was reexamined recently by 
Bahcall, Steinhardt, & Schlegel (2004). Other recent attempts include measurements of 
absorption hue spectra in alkah-like atoms (Potekhin & Varshalovich 1994, Varshalovich 
& Potekhin 1994, Murphy et al. 2001c, Fiorenzano, Vladilo, & Bonifacio 2003). While 
many observations have led to interesting limits on the temporal variation of a (for a recent 
review, see Uzan 2003), only the MM method has led to a quantitatively positive result. 
Murphy et al. (2003a) deduce that 6a/a = (—0.54 ± 0.12) x 10~^ over a redshift range of 
0.5 < z < 3, where 6a is defined as the deviation from present value. The implications 
of this deduced variation in a at around a 5a significance are phenomenal, and several 
cosmological models to explain its origin have been proposed (see e.g., Beckenstein 1982, 
Sandvik, Barrow, & Magueijo 2002, Olive & Pospelov 2002, Wetterich 2003, Anchordoqui 
& Goldberg 2003, Copeland, Nunes, & Pospelov 2004, Lee, Lee, & Ng 2003; Byrne & Kolda 
2004). The caveat of implementing such a precise method, however, is its possible sensitivity 
to unknown systematic errors. 

Some of the excitement concerning the evolution of the fine structure constant has 
been focused on finding alternative explanations of the observed line structures or other 
systematic errors. Chand et al. (2004) and Srianand et al. (2004) probed the sensitivity 
of the MM method with respect to synthetic signal alterations and found that the MM 
method may break down in well blended, multi-cloud systems. They also applied the MM 
method independently and found 6a/a = (—0.06 ± 0.06) x 10^^. Another group (Quast, 
Reimers, & Levshakov 2004) has also recently applied the MM method utilizing exceptionally 
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high-resolution QSO absorption-line spectra. Their results are also consistent with a null 
hypothesis regarding the fine structure evolution. Other potential systematic errors in the 
MM method have been elicited by others. They involve the cloud velocity structure and line 
blending (Bahcall et al. 2004) or cloud inhomogeneity and spectrographic inconsistencies 
(Levshakov 2003). 

A number of sources of possible systematic error in this method have been well docu- 
mented (Murphy et al. 2001b and 2003b; see also Bahcall et al. 2004). Here, we will focus 
on one of these possible systematic errors for which there is recent evidence of a new inter- 
pretation, namely the isotopic abundances of Mg assumed in the analysis. In this paper, we 
expand on our earlier studies (Ashenfcltcr, Mathews & Olive 2004) of possible systematic 
effects from the chemical evolution of magnesium isotopes within DLA quasar absorption- 
line systems. All of the results quoted above are based on the assumption that the isotopic 
abundance ratios of Mg are Solar in these systems. Based upon galactic chemical evolution 
studies previously available (Timmes, Woosley & Weaver 1995), one could argue that the 
ratio of ^^'^^Mg/^^Mg is expected to decrease at low metallicity. In this case, if it is assumed 
that only ^^Mg is present in the absorbers, the Murphy et al. (2003a) result becomes signif- 
icantly stronger 6a/a = (—0.86 ±0.10) x 10^^ (assuming also that only ^^Si is present) and 
the Chand et al. (2004) limit becomes a detection 6a/ a = (—0.36 ± 0.06) x 10"^. Hence, it 
is possible that the detections of time- varying a are even more significant than the quoted 
confidence limits. However, we show that it is also plausible that the ^^'^^Mg/^^Mg ratio 
was in fact sufficiently higher at low metallicity to account for the apparent variation in a as 
seen in the so-called "low redshift" (0.5 < z < 1.8) data. Thus, the MM method of analysis 
may provide important new insights into the chemical evolution of damped Lyman-a quasar 
absorption-line systems rather than conclusive evidence for a time-varying fine-structure 
constant. 

We begin the present discussion by briefiy reviewing the current observational limits 
on the variations of the fine structure constant. In section 3, we discuss the theory and 
observations of the Mg isotopes. The sensitivity of {6a /a) to the Mg isotopic abundances is 
explained in section 4. In section 5, we describe a simple chemical evolution model, which we 
then utilize to address the question of the history of the Mg isotopes and other elements as 
a function of metallicity. Results of this study and the sensitivity to the model are discussed 
in section 6. Our summary and conclusions are given in section 7. 
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2. Limits on the fine-structure evolution 

There are a number of important astrophysical and terrestrial constraints on the fine- 
structure constant that must be respected. The most primordial of the limits comes from big 
bang nucleosynthesis (Kolb, Perry, & Walker 1986, Malaney & Mathews 1993; Scherrer & 
Spergel 1993, Campbell & Olive 1995, Bergstrom, Iguri, & Rubinstein 1999, NoUett & Lopez 
2002), which tests for variations back to a redshift as high as ~ 10^°. However, the limit 
attained is rather weak {5a/ a) < 10"^. If one assumes that variations in the fine structure 
constant are accompanied by variations in other gauge and Yukawa couplings (Campbell & 
Olive 1995, Langacker, Segre, & Strassler 2002, Dent & Fairbairn 2002, Calmet & Fritzsch 
2003, Damour, Piazza & Veneziano 2002), this limit can be strengthened by about two orders 
of magnitude (Campbell & Olive 1995, Ichikawa & Kawasaki 2002). There is also a slightly 
weaker constraint {6a /a) < few xlO~^ from the WMAP cosmic microwave background power 
spectrum corresponding to the epoch of photon last scattering at a redshift of 2; ~ 1100 
(Rochaet al. 2003). 

At smaller look-back times (lower redshifts), there are significantly stronger limits. Me- 
teoritic data on the radioactive /9-decay of ^^^Re was used to place an upper bound of 
Sa/a < 10"^ (Olive et al. 2002, 2004, Fujii & Iwamoto 2004) and is applicable to a redshift 
of 2; < 0.45. This limit improves by a factor of ~ 25 when variations in a are assumed 
to be coupled with variations in other physical constants. The strongest limit is based on 
the ^^''Sm resonant neutron-capture cross-section operating in the Oklo natural fission reac- 
tor (Shlyakhter 1976, Damour & Dyson 1996, Fujii et al. 2000) about 2 billion years ago 
{z ~ 0.15). The limit on a is 5a;/a; < 10^^ and can be improved by 2-3 orders of magnitude 
when coupled with variations in other constants (Olive et al. 2002). Finally, atomic clocks 
provide very stringent constraints on the present-day rate of change of a. By comparing hy- 
perfine transitions in ^^Rb and ^^^Cs using atomic fountains over a period of 5 years, Marion 
et al. (2003) were able to derive the limit a/a < 1.5 x 10"^^ yr~^. Combined with data from 
Hg+ and H, Fischer et al. (2003) obtain a shghtly stronger bound 5a/a = (-0.9±4.2) x 10"^^ 
over a 3-4 year period. For a comprehensive review, see Uzan (2003). 

Interestingly enough, none of the limits above cover the redshift range corresponding to 
the quasar absorption-line DLAs that yield the recent evidence of a variation in a. Hence, 
the door remains open for creative model building. Nevertheless, the strength of these 
limits, particularly the Oklo and ^^^Re bounds, may indicate that something other than a 
time-varying a may be responsible for the effects seen in the many-multiplet method. 
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3. Mg: Theory and Observations 

3.1. Expected trends of Mg production 

The MM method is sensitive to the Mg isotope abundances for data in the 0.5 < 2; < 
1.8 redshift range for which the apparent variation in a is most pronounced (Murphy et 
al. 2003b). Hence, it is imperative to carefully scrutinize the chemical-evolution history that 

one expects to occur within the QSO absorption systems. The DLAs that were studied 
to obtain the apparent variation in the fine-structure constant are likely to be galaxies in 
various stages of evolution. Fenner, Prochaska & Gibson (2004) conclude that the sources of 
DLAs span a wide range of galaxy morphological types and sizes, from dwarf irregulars to 
giant ellipticals. Hence, chemical evolution models that fall within these broad morphologies 
can be used to explore the expected trends in the abundance distributions. 

Mg is produced in both Type 1 and Type 11 supernovae. In Type 11 supernovae, it 
is produced in the carbon and neon burning shells with an abundance somewhat less than 
10 percent of the oxygen abundance produced in massive stars (Woosley & Weaver 1995; 
henceforth WW95). ^^'^^Mg are produced primarily in the outer carbon layer by the reactions 
2^Ne(Q;,n)25Mg and 2^Mg(n,7)26Mg. The Solar-metaUicity models of WW95 eject Solar values 
with Mg isotopes ratios reasonably close to the terrestrially observed value of ^^Mg:^^Mg:^^Mg 
= 79:10:11 (Rosman & Taylor 1998). More massive stars tend to be slightly enhanced in 
the heavy isotopes (eg., the WW95 25 Mq model gives a ratio of 65:15:20). Furthermore, 
the abundance of ^^'^^Mg scales linearly with metallicity in the carbon shell. As a result, 
it would be expected that the ejecta from the first generation of Type-11 supernovae would 
contribute almost no ^^'^^Mg. Some of the Solar abundance of Mg is produced in Type I 
supernovae but with the [Mg/Fe] far below Solar. For example, the models of Thielemann, 
Nomoto, & Yokoi (1986) give [Mg/Fe] = -1.2. Due to the absence of free neutrons, essentially 
no ^^'^^Mg is produced in Type I supernovae. Thus, in previous chemical evolution models 
the heavy Mg isotopes could only be efficiently produced at late times in high metaUicity 
Type H supernovae. 

Moreover, of critical importance in the present study is that ^^'^^Mg can also be produced 
in intermediate-mass (IM) asymptotic-giant-branch (AGB) stars. It has been recently noted 
(Karakas & Lattanzio 2003; Siess, Livio, & Lattanzio 2002; Forestini & Charbonnel 1996) 
that IM stars of low metallicity can in fact be efficient producers of the heavy Mg isotopes 
during the thermal-pulsing AGB phase. 

It should be noted that the Mg yields of WW95 near 20 Mq are substantially lower 
than the predictions of other stellar models as is described in Argast, Samland, Thielemann, 
& Gerhard (2002). While this underestimate only affects a narrow mass range in the Type 
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II yields, it may have an impact on the total Mg yield in IMF weighted chemical evolution 
models. If the '^'^Mg production from the Type-II models were corrected upward with respect 
to the WW95 yields, the expected heavy isotopic ratio at low metallicity would decrease. 
25,26jYfg scaigg "with the metallicity, so it would not be affected by this uncertainty. The 
primary consequence of the "^^Mg underestimate in this mass range would be to slightly 
enhance the impact that the Mg yields from AGB stars on the Mg isotopic ratio. However, an 
enhancement of the Type-II Mg yields could be compensated by an appropriate adjustment 
to the IMF as described below in Section 5. 



3.2. Mg production in AGB stars 

Heavy magnesium isotopes arc synthesized via two mechanisms both of which arc par- 
ticularly robust in 2.5 - 6 Mq stars with low metallicity. Such low-metallicity stars are 
precisely the kinds of objects which ought to produce the abundances observed in damped 
Lyman-cc QSO absorption-line systems at high redshift for the reasons already described. 

One process particularly effective in low-metallicity AGB stars (Boothroyd, Sackmann 
& Wasserburg 1995) is that of hot-bottom burning (HBB). During the AGB phase, stars 
develop an extended outer convective envelope and material in the hydrogen envelope is 
mixed downward to regions of high temperature at the base. Of particular interest is that the 
base of the envelope is more compact and at higher temperature in low-metallicity stars than 
in stars of Solar composition. This can be traced to the decreased opacity of these objects. 
Because these stars become sufficiently hot (T > 7 x 10^ K), proton capture processes in the 
Mg-Al cycle become effective. Proton capture on ^"^Mg then leads to the production of ^^Mg 
(from the decay of ^^Al) and to ^^Al (which decays to ^^Mg). The relevant thermonuclear 
burning reactions have a strong sensitivity to temperature. Hence, IM progenitor stars of 
low metallicity can contribute abundant products from high-temperature burning to the 
interstellar medium. Moreover, these stars should also have a shorter timescale than their 
higher-metaUicity counterparts, though this effect is not large as 5-6 Mq stars have relatively 
short lifetimes at any metallicity. 

A second contributing process occurs deeper in the star during thermal pulses of the 
helium-burning shell. The helium shell experiences periodic thermonuclear runaways when 
the ignition of the triple-alpha reaction occurs under electron-degenerate conditions. Due to 
electron degeneracy, the star is unable to expand and cool. Hence, the temperature rapidly 
rises until the onset of convection to transport the energy away. During these thermal pulses, 
^^Ne is produced by a-capture on ^^N, which itself is left over from the CNO cycle. Heavy 
magnesium isotopes are then produced via the ^^Ne(Q!,n)^^Mg and ^^Ne(Q;, 7)^^Mg reactions. 
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Note that this process is important for stars with M Ri 3 M© (Karakas & Lattanzio 2003), 
while HBB is more important for M ^ 4 M©. 

Several groups give credence to the assertion that AGB stars produce heavy Mg isotopes. 
Siess, Livio, & Lattanzio (2002) describe the ^^'^^Mg production that occurs during the 3rd 
dredge-up. A key point is that even though seed material is less plentiful in low-metallicity 

stars, the reactions are very temperature sensitive. Hence, the increased temperature in 
the interior of low-metallicity stars more than compensates for the depleted seed material, 
leading to significant production of the heavy Mg isotopes. Moreover, seed material produced 
and mixed during the first two dredge-up episodes will also be more efficiently produced due 
to the heightened temperature. It has even been argued that these processes may also be 
net destroyers of ^^Mg (Forestini & Charbonnel 1996, Karakas & Lattanzio 2003) due to the 
extreme temperatures attained. 

To illustrate this phenomenon regarding the Mg evolution in AGB stars, Denissenkov & 
Herwig (2003) modeled a typical low-metallicity {Z = 0.0001) AGB star of 5 M© and found 
that it was capable of taking an initial ratio of ^^Mg:^^Mg:^^Mg = 90:5:5 to the extreme 
ratio of (13:71:16)! It should be noted that a Mg isotopic composition this enriched in ^^Mg 
has never been detected. Nevertheless, their model estabhshes an upper limit on the effect 
of HBB on material dredged up to the surface. The main impetus of their study was to 
qualitatively illustrate that neither low seed material nor recently updated reaction rates 
can prevent the efficient processing of Mg in AGB stars. 



3.3. Observations of Mg abundances 

The data on Mg abundances in low-metaUicity stars exhibits considerable dispersion. 
This dispersion is in excess of a factor of 3 for a fixed value of [Fe/H]. While such dispersion 
could be a symptom of systematic error when data from several samples are combined, the 
observed intrinsic dispersion in abundances of low-metallicity stars is generally interpreted 
(e.g. Ishimaru & Wanajo 1999) as evidence for effects of local stochastic star-formation 
events. (We note, however, that the recent data of Cayrel et al. (2004) show considerably 
less dispersion in [Mg/Fe] and other [a/Fe] ratios in halo stars than previous results. Hence, 
the need for highly inhomogeneous halo evolution is correspondingly reduced.) 

Gay & Lambert (2000) determined the Mg isotopic ratios in 20 stars in the metallicity 
range —1.8 < [Fe/H] < 0.0 with the aim of testing theoretical predictions (e.g. Timmes et 
al. 1995). These data exhibit a large dispersion in the ^^'^^Mg/^^Mg isotope ratio. Hence, 
the case for a low ^^'^^Mg/^^Mg ratio at low [Fe/H] was, perhaps, not yet unambiguously 
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established, although their results indicated that ^^'^^Mg/^^Mg appears to decrease at low 
metallicity for normal stars. Nevertheless, based upon the large dispersion in these data, 
one could neither make the case for a high nor low ratio of heavy Mg isotopes. However, 
indirect evidence for Mg/Al-cyclc clement production in low-metallicity stars is suggested 
by the observations of Denisscnkov ct al. (1998), who found substantial Al enhancements in 
globular-cluster giants at the expense of Mg abundances. 

Although many of the stars studied by Gay & Lambert (2000) were found to have Mg 
isotopic abundance ratios somewhat higher than predicted, even the "peculiar" stars which 
show enrichments in ^■''•^^Mg do not have abundance ratios substantially above Solar. More- 
over, it should be noted that Timmes et al. (1995) point out that their model underestimates 
the ^^'^^Mg by around a factor of 2 at a metallicity of [Fe/H] = — 1, and that this discrepancy 
may indicate that another source of magnesium isotopes is in operation at low metallicity. 

Based on the galactic chemical evolution models available at that time (e.g. Timmes 
et al. 1995), the adoption of Solar isotopic Mg ratios by Murphy et al. (2003a) in the MM 
method would appear to have been safe and conservative. Previous models seemed to indicate 
that the ^^Mg:^^'^^Mg ratio was higher than 79:21 in the past, and it was shown that a 
composition more rich in ^^Mg only strengthens the case for the evolution of a. However, 
as we discuss in more detail below, increasing the abundances of the heavier Mg isotopes 
would yield a larger multiplet splitting from isotope shift effects. This would imply an 
apparently higher value for a (i.e. closer to the present value). Indeed, raising the heavy 
isotope concentration to ^"^Mgi^^'^^Mg = 63:37 could remove the significance of the signal for 
a time- varying the fine-structure constant. 

In support of this possibility, a new study of Mg isotopic abundances in stars in the 
globular cluster NGC 6752 (Yong et al. 2003) affects the case for assuming a Solar isotopic 
composition. This study looked at 20 bright red giants with an assigned metallicity of [Fe/H] 
= -1.62. Since globular clusters may be the remnants of early galactic chemical evolution, 
their abundances might be representative of the abundances to be found in the DLAs of 
interest in the present study. The observations of Yong et al. (2003) show a considerable 
spread in the Mg isotopic ratios which range from ^^Mg:^^Mg:^^Mg = 84:8:8 (shghtly poor 
in the heavies) to 53:9:39 (greatly enriched in ^^Mg). Of the 20 stars observed, 15 of them 
show ^^Mg fractions of 78% or less (i.e. enriched in heavy isotopes relative to Solar), and 7 
of them show fractions of 70% or less with 4 of them in the range 53 - 67 %. This latter 
range, if representative of the abundances in the DLA systems, is sufficiently low to have 
substantially affected any determination of a. A previous study (Shetrone 1996) also found 
unusually high abundances of the heavy Mg isotopes in giant stars in the globular-cluster 
M13. For this system, ^^Mg:^^'^^Mg as low as 50:50 and even 44:56 was found. Such values 
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exceed the necessary condition to account for the apparent variation in the fine structure 
constant in the low redshift sample of Murphy et al. (2003). 

Other chemical evolution models of low-metallicity systems also suggest that there may 
be a missing contributor to the cosmic Mg abundance. Standard chemical evolution models 
were compared to observations of DLAs by Fenner et al. (2003; 2004). These studies indicate 
that supplemental sources of Mg beyond massive stars are needed. They concluded that, 
although Mg is dominated in the present-day ISM by SN ejecta, SNe probably played a lesser 
role to that of IM stars in earlier epochs. 

It should be emphasized that the a derived from the MM method is sensitive to the 
isotopic ratio of Mg and not the total abundance. The missing Mg from chemical evolution 
models and abundance measurements points to another source, one which played a dominant 
role in the earliest epochs of star formation. In our models, this source is derived from AGB 
stars. 



4. Sensitivity of the MM method to Isotopic Shifts of Mg 

Before we describe the specifics of the chemical evolution model, it is necessary to illus- 
trate how the heavy isotopes affect the determination of a from the MM method. Normally, 
the MM method would compare the line shifts of the species particularly sensitive to any 
real change in a to one with a comparatively minor sensitivity (referred to as an anchor). 
For the low redshift sample, which consists of 74 out of 128 measured absorption systems 
[and the most significant evidence for variation in the Murphy et al. (2003a) data as well 
as all of the data in the recent analysis of Chand et al. (2004)], Fe fines are compared to 
Mg lines (which serve as the anchor). For our purposes, the sensitivity of the MM method 
to the fine structure constant can be roughly approximated by subtracting the wavenumber 
shift of Fe from Mg 

/^UJMg - AtUFe = ?e//^ = ge//2 , (1) 

a 

where Aa;i is the difference between the observed wavenumber of line i and the laboratory 
value. The quantity X relates to the change m. a, X = (a^/ao)^ — 1 f« 2^, where the 
approximation is adequate for the small variations of interest here. While the true method 
simultaneously fits ratios of Fe lines to Mg lines, the effects of isotope shifts can be estimated 
by taking the average trend as is described below. The quantity qeff is the difference between 
the average value of q for Mg and that for Fe. From the q values given in Murphy et 
al. (2003a), we obtain q^ff ~ 1280 ± 150 from an average of the 74 low- 2; absorption systems. 
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The same wavenumber shifts can be accomplished from the isotopic shift (IS) of Mg 
alone (since the IS for Fe is small) 

AuJMg - Aco-Fe = ^A' - ■ (2) 

The IS depends sensitively on the field shifts, specific mass shifts, and normal mass shifts 
that in turn depend upon the nuclear charge configuration. Hence, it must be determined 
experimentally. Accordingly, Berengut, Dzuba & Flambaum (2003) have experimentally 
deduced the coefficients of these shifts. They also provided a relationship between isotopic 
abundances and frequency shifts, 

UJA' -OJa^ -^[{Kms + Kms){^, - ^) + F5 < P >^''^] , (3) 

where A' is the mean atomic mass (enhanced with the heavier isotopes) and A is the Solar 
atomic mass number (for Mg, A = 24.32). The quantity ksms is the specific mass shift 
coefficient, while knms is the normal mass shift. F is the field shift, while 5 < >'4',a 
is the difference in the mean square radius, and c is the speed of light. Incorporating the 
coefficients given in Berengut et al. (2003), we can determine the corresponding effective 
a-variation for the A' corresponding to a specific isotopic composition of Mg. Recall that 
increasing the abundances of the heavier Mg isotopes yields a larger value for the deduced a, 
and a ratio of ^^Mg:^^'^^Mg = 63:37 is sufficient to obviate the need for a time-varying fine 
structure constant. Note that the IS also depends on the relative ratio of ^^Mg versus ^^Mg. 
The heavy isotopic ratio of (25Mg+2^Mg)/2^Mg = 37/63 = 0.58 quoted above still assumes 
a Solar ratio of ^^Mg/^^Mg. 

The change in the deduced values of 5a/ a relative to those obtained when assuming 
a Solar isotopic composition (i.e. the correction to 5a/ a due to non-Solar values of the 
Mg isotopic abundances) is shown in Fig. 1. Therefore, for Solar isotopic composition 
[(^^Mg+^^Mg)/^'^Mg = 0.27], the value of a is the same as that obtained by Murphy et 
al. (2003a) and Chand et al. (2004). When only ^^Mg is assumed to be present, there is a 
shift if about —0.67 x 10~^, thereby amplifying the signal they would have reported. Note 
that this shift is somewhat stronger than the shift quoted in Murphy et al. (2003a) as we 
are using the more recent results of Berengut, Dzuba & Flambaum (2003) and Berengut et 
al. (2003). 

In contrast, when the ratio of heavies to ^"^Mg is larger, the deduced 5a is smaller and 
can effectively cancel the effect obtained by Murphy et al. (2003a). Using our adopted mean 
value of qeffi we find that the Mg isotopic ratio required to cancel the apparent result of 
Murphy et al. (2003a) is {^^Mg-i-'^^Mg) /'^'^U.g = 0.62. Figure 1 also shows the dependence of 
5a on the ratio of ^^Mg:^^Mg in addition to the total isotopic ratio. If the ratio ^^Mgr^^Mg 
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Fig. 1. — Effects of the Mg isotope shift. Plotted is the change in deduced values of 5a/ a 

relative to values obtained from an assumed Solar Mg isotopic composition. The change 
in the deduced value is plotted as a function of the (^^Mg+^^Mg) /^"^Mg isotopic abundance 
ratio. The solid curve is the value deduced for a Solar ratio of ^^Mg "with respect to ^^Mg. 
The dotted (dashed) curve corresponds to ^^Mg:^^Mg = 2:1 (3:1). The dot-dashed curve 
shows a sub-Solar ratio of ^^Mg:^^Mg = 8:10. The horizontal lines mark the value of 5a/ a 
needed to compensate for the apparent shift in a determined by Murphy et al. (2003a) 
(dashed) and Chand et al. (2004) (thin sohd). 

were 2:1 or 3:1, the necessary total isotopic ratio, (^^Mg+^^Mg)/^''Mg, would be only 0.54 
and 0.50, respectively. A sub-Solar ratio of ^^Mg/^^Mg = 8:10 would require the total isotopic 
ratio to be 0.68. Fig. 1 also illustrates that the ratio of ^^Mg:^^Mg has a pronounced effect 
on the deduced value of 5a/ a. The deduced vahic of 5a/ a increases by ~ 1.5 x lO^*' even 
for a Solar ratio of (^^Mg-|-^^Mg)/^''Mg = 0.27 just by varying the ratio of ^^Mg:^^Mg from 
1:1 to 3:1. 

The change in the deduced value of 5a/ a is also sensitive to the value of Qeff- While 
we have approximated the MM method with an cfTcctivc average from the low-z absorption 
systems, this value is somewhat uncertain. In Fig. 2, we illustrate the sensitivity of the 
deduced value for 5a/ a to ge//- Adopting a Solar ratio for ^^Mg:^^Mg and taking the la- 
range of qeffi we obtain a mean value and uncertainty in the Mg isotopic ratio needed to 
compensate the apparent deviation in the fine structure constant reported in Murphy et 
al. (2003a) of (25Mg+26Mg)/24Mg = 0.62 ± 0.05. Similarly, we infer (25Mg+26Mg)/24Mg = 
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Fig. 2. — The shift in deduced values of 5a/ a relative to values deduced from an assumed 
Solar Mg isotopic composition resulting from different values of (/,//• The solid line is the 
shift using the mean value of Qeff = 1280. The dotted lines sho'w the effect of the la limits 
on Qeffi ±150, "where the uncertainty is dominated by the uncertainty of the qp^ given in 
Murphy et al. (2003a). The steeper (flatter) slope corresponds to ^e// = 1130 (1430). 

0.30 ± 0.01 for the Chand et al. (2004) data. This uncertainty in gfe// translates into an 
uncertainty in 5a/ a of ±0.14 x 10~^ near a shift in 5a/ a near 5.4 x 10~^. 

While 5a/ a is sensitive to the Mg isotopic composition in the low redshift sample, the 
systems at 2; > 1.8 primarily depend upon Si. In Berengut et al. (2003), it "was determined 
that the Si// lines are relatively insensitive to a heavy isotopic shift, although heavy Si 
isotopes may be produced in an analogous fashion to the Mg isotopes in AGB stars. Fully 
accounting for the 5a/a in the sample with z > 1.8 via isotopic abundances would require 
(Berengut et al. 2003) a conspiracy of several isotopic abundances. 

5. DLA chemical evolution with early AGB enrichment 

In order to explore the Mg isotopic evolution in DLAs, it is important to build a plau- 
sible chemical evolution model that is constrained by observations. Unfortunately, detailed 
abundances for most of the systems observed in the Murphy et al. (2003a) have not been 
reported. Nevertheless, we can assume that their sample of DLA absorption systems have 
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similar characteristics to those of measured DLA systems within the same redshift range. 
Systems in the redshift range of interest here (0.5 < z < 1.8) typically span a broad range 
of metallicities with the mean value around [Fe/H] = -1.1 (Pettini 1999, Cen et al. 2003). 
Pettini (1999) also demonstrated that there is no discernible evolution of metallicity in the 
redshift range oi z = 0.4 — 3.5. This lack of metallicity evolution in DLAs is further sup- 
ported by Prochaska & Wolfe (2000, 2002). DLAs are distinguished by their high neutral 
hydrogen column densities; however, if a system were to produce large amounts of metals, it 
would also tend to ionize the neutral hydrogen gas. Alternatively, systems with high column 
densities in neutral hydrogen should also have reached the conditions for star formation. In 
fact, no observed DLA has a metallicity lower than [Fe/H] < —3 (Cen et al. 2003). Hence, 
we confine our chemical evolution models to the broad range of metallicities observed in 
DLA systems. 

In this section, we describe a chemical evolution model that can account for the observed 
variation in the fine structure constant. We utilize recent Mg isotopic yields from AGB 
stars. Wc then compare the model abundances to observed abundances in DLA systems as 
a constraint on this interpretation of the MM results. 

For our purposes, a simple recalculation of the model of Timmes et al. (1995) with 
and without the enhanced contribution from AGB stars is sufficient. This allows us to 
make a direct comparison with the conclusions of the previous authors. We utilize the same 
metallicity-dcpcndcnt yields of WW95 for Type II supernovae, using the same choices for 
explosion energies as in Timmes et al. (1995), [i.e. model A of WW95 for m = 11 — 25 Mq, 
and model B for m = 30, 35, and 40 Mq, which was taken as the upper limit of the IMF.] 
We adopt the metalhcity-dcpcndent stellar lifetimes and the yields for the stellar mass range 
from 6 — 7MQ given by Portinari, Chiosi, & Bressan (1998). We supplement the IM yields 
of Marigo (2001) for stellar masses up to 5M0 (using the a — 2 mixing length parameter 
model) with the recent AGB yields of Karakas & Lattanzio (2003) where the abundances of 
all three Mg isotopes are included. 

Although the yields for massive stars are available from zero to Solar metallicity, we 
need to extrapolate the yields for IM stars below the lowest adopted metaUicity model of 
Z — 0.004. We keep the yields constant when below the lowest metallicity model. This 
approach is adequate for primary elements, but needs to be justified for elements such as N 
and Mg, which are known to derive from secondary sources at high metallicity. Although 
N is largely a secondary element at high metallicity, it is known (Pagel 1997) to be largely 
primary at low metallicity (Z^O.004). Hence, constant nitrogen for Z^O.004) seems justified. 
Likewise, Mg yields can be held constant, since it roughly scales with N. Obviously, these 
assumptions leads to some uncertainty and needs to be checked with future calculations of 
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IM stars at lower metallicity. 

In addition to extrapolationg to low metallicity, it was also necessary to interpolate in 
the mass range of 7 — 12 Mq between the IM models and WW95 models. This leads to an 
uncertainty in the yields for that mass range. A simple linear interpolation between the IM 
and WW95 models would be dominated by the larger yields of the Type II SNe. Also, stars 
in that mass range at some point no longer experience the thermally pulsing AGB phase. 
Therefore, we adopt a more conservative approach whereby we interpolate the yields for both 
models to zero for 9 Mq stars in order to distribute the uncertainties in the yields between 
both sets of yield models. We also adopt the Solar abundances of Grevesse & Sauval (1998). 

The model of Timmes ct al. (1995) is based upon an exponential infall of primordial gas 
over a timescale of 4 Gyr, and a Schmidt law for the star formation rate. For our purposes, 
the efficiency of star formation can be modified to best account for the observed abundances 
in DLAs. This is particularly true given that they are most likely the progenitors of a wide 
range of galaxy morphologies. 

We assume instantaneous mixing with no outflow. Hence, the evolution of the mass 
fraction Xi of isotope i can be written. 



where, B{t) (in units of Mq Gyr"^) is the stellar birthrate function at time t, ^(m) is the 
initial mass function (IMF), Xf is the mass fraction of isotope i ejected by various sources 
specified by superscripts {x — S for normal evolution, x — la for supernovae Type la ejecta, 
and X — AGB for the supplement ejecta of Karakas & Lattanzio(2003)). The lifetime of 
a progenitor star of mass m is denoted T(m), while mco is the mass of the carbon-oxygen 
white-dwarf progenitor of the Type la supernovae, with Ria as the rate of Type la supernovae. 
The quantity Mgas is the gas mass at time t. Mg^i represents the galactic infall rate of isotope 
i (presumed to be primordial material). The third term on the right hand side of Eq. 4 
refers to the amount of isotope i incorporated into new stars. The yields given by Karakas 
& Lattanzio (2003) are for AGB stars up to 6.5 Mq. Even at m = 6.5 Mq the production 
of heavy Mg isotopes is still increasing with mass. Therefore, we extrapolate the Karakas & 
Lattanzio (2003) yields to 7 Mq. Beyond that, the effective Mg yield is linearly interpolated 
to zero at 8 Mq because such stars do not ignite He shell burning under electron-degenerate 





= mcoXl^Ria + Mg^i - B{t)Xi{t) 




(4) 
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conditions. Neither do they go through an AGB phase. Although our results are somewhat 
sensitive to this extrapolation, a change in our result due to a different extrapolation could 
be compensated by a change in the AGB enhanced IMF. 

While the stellar evolution models of Karakas & Lattanzio (2003) show robust AGB 
evolution beyond 5 M0, the adopted CNO yields of Padova stellar models of Marigo (2001) 
and Portinari et al. (1998) do not exhibit AGB evolution beyond 5 M0. Until self- consistent 
yield models of AGB stars are developed that follow both Mg and CNO evolution, we are 
limited to supplementing Mg evolution in the manner previously described. Consequently, 
this approach may affect the direct coupling of the Mg to the associated CNO abundances. 

In our model, the cosmic Type la supernova rate is given by the formulation of Kobayashi, 
Tsujimoto, & Nomoto (2001). This rate incorporates a minimum metallicity condition be- 
fore the operation of Type la supernovae of [Fe/H] > —1.1. This is adopted as a necessary 
condition for the white dwarf progenitor to accrete effectively from the binary companion. 
This condition was hypothesized by Nomoto et al. (2003) based upon the fact that the stellar 
wind velocity from the binary companion is slower at low metallicity and also that the mass 
range of white-dwarf progenitors is somewhat limited for the short timescales corresponding 
to low metallicity. 

We propose a modest early enhancement of the IMF along the lines of Fields et al. (1997; 
2001). We present some theoretical and observational evidence for this enhancement below. 
Even though the early IMF is enhanced in the IM range, the usual normalization is still 
adopted 



^{m)dm — 1 (5) 

We then write the star formation rate B"^ as 

B{t)^{m) = Bi{t)^i{m) + B2it)^2im) (6) 



= Sim"^-^' + {B2/m)exp -{\og{m/mc) )/{2a^) 

The ^i(m) IMF accounts for a standard Salpeter distribution of stellar masses. ^2(^77.) is 
an additional log-normal component of stars peaked at rric Mq. The dimensionless width, 
a, allows the mass distribution to extend across the entire IM range. For the model of 
Ashenfeltcr, Mathews, & Olive (2004), nic = 5 Mq and a = 0.07 was adopted. In the 
present work, a range for these parameters is adopted consistent with the constraints deduced 
in Adams & LaughUn (1996). 

For the normal stellar component, we parametrize the time dependence of the stellar 
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birthrate function Bi{t) as 



2 



= (1.0 - e-'/^')esFMtot{t) MjMtot{t) 



(7) 



while for the IM component we similarly adopt 



B2{t) = A2e-"'''esFMtot{t)[ 



2 



(8) 



Mtot{t) 



This model thus contains an early burst of IM stars peaked at rric, with a spread of 
masses governed by the dimensionless width a. The burst is exponentially suppressed on a 
time scale of T2. The Bi component describes the standard quiescent star formation with 
a smooth transition from the burst. We note that the IM component obviously dominates 
the mass recycle rate early on (depending on the values for the free parameters) , but quickly 
gets diluted as the normal Salpeter component evolves. The coefficient A2 in B2{t) was 
adjusted to produce sufficient Mg isotope enhancement. Finally, esF is the coefficient of star 
formation efficiency. For reference, Timmes et al. (1995) determined that the best fit for the 
Solar neighborhood was esF = 2.8. 

It should be noted that this coefficient gets normalized with the IMF, so that it only 
affects the abundance ratio when the log-normal enhancement is competing with the normal 
power law component of the IMF. The normalized IMF is shown at four different times in 
Fig. 3. The sensitivity of our results to each of the parameter choices in Eqs. (6) - (8) will 
be discussed in section 6 below. 

In addition to the overall star-formation efficiency, egFi the model consists of two pa- 
rameter categories: the parameters that determine the total amount of IM stars above the 
standard IMF (A2, Ti, T2) and those that determine the distribution within the enhancement 
(mc, cr). The exponential decay of the enhancement of IM stars ensures that a smooth tran- 
sition into the standard IMF occurs. The combination of A2 and T2 determines the ratio of 
the two IMFs while they are in competition. 



To some extent, our adopted increase in the IMF for IM stars may simply be thought of 
as a compensation for the uncertainties in the nucleosynthesis yields of AGB stars. Indeed, 
a very small change in burning temperature or a more efficient dredge-up could easily ac- 
complish the required enhanced ^^'^^Mg yields. Nevertheless, there are also reasons to expect 
that the early IMF would be considerably different than the present day IMF. Hence, we 
adopt an enhanced AGB-star IMF as a reasonable parameterization. 



5.1. Arguments for an AGB-enhanced IMF at low metallicity 
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Fig. 3. — Evolution of our adopted AGB-enhanced IMF. The log-normal component is 
centered at 5 Mq with a = 0.07 as in our model 1. This enhancement, however, quickly 
evolves to a standard Salpeter IMF. The values of ri and T2 are 0.5 and 0.2 Gyr, respectively. 

Arguments for an enhanced IMF for IM stars are as follows: with fewer metals initially, 
cooling in the dense protostellar cloud is predominantly done by atomic hydrogen, which 
is less efficient. Therefore, a more massive cloud is required to cool and gravitationally 
collapse into stars. The lower limit of the IMF depends sensitively on the length scale of 
density fluctuations and turbulent fragmentation. These density fluctuations in the present- 
day interstellar medium (ISM) are largely the result of previous star formation, which was 
obviously absent in the first star formation epoch. Moreover, these primordial clouds are at 
higher temperature in the low-metallicity ISM. Therefore, they further inhibit the smaller 
length scales of the density fiuctuations. Uehara et al. (1996) set the lower mass fimit of the 
extremely metal poor stars at IAMq, which consequently implies that all metal- free stars 
have since evolved into white dwarfs or other remnants. Although a 1.4 Mq cut-off does not 
necessarily imply a higher relative abundance of IM with respect to massive stars, arguments 
of a more steeply sloped early IMF have also been proposed which would suppress massive 
stars relative to the IM stars of interest here. For example, Padoan & Nordlund (2002) 
suggest a steeper IMF slope as due to weaker magnetic fields in molecular coulds at high 
redshift. 

A further constraint is provided by the models of Yoshii & Saio (1986), who concluded, 
based on an opacity-limited fragmentation model, that a metal-free IMF should be peaked at 
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4 — IOMq. It should then dedine steeply on both sides of the peak. Furthermore, competing 
models of a metal-free IMF have the peak at much higher mass (Hernandez & Ferrara 2001) 
and can even allow for very massive objects (VMOs) to form (Susa, Uehara, & Nishi 1996). 

Further potential constraints on an enhanced IMF for IM stars exist from observations 
of stars in the Galactic halo, although we emphasize that details of the local Galactic envi- 
ronment do not necessarily apply to the Lyman-a systems of interest to the present study. 
Even so, Adams & Laughlin (1996) and Fields et al. (1997) (see also Ryu, Olive & Silk 
1990) delved into the form of the IMF that could match the MACHO data of gravitational 
microlenscs in the Galactic halo. They independently determined that the early IMF is 
tightly constrained to be in the mass range of 1 < Mq < 8. These IM stars would evolve 
into white dwarfs and populate the halo if they had formed early enough to evolve and cool. 
Specifically, they determined that the maximum of the log-normal distribution is around 
2 — 3Mq and is sharply peaked. 

Regarding those earlier models that produce a very strong enhancement of IM stars, 
it has since been argued that those models overproduce current constraints on halo white 
dwarfs (Gibson & Mould 1997; Flynn, Holopainen & Holmberg 2003; Garcia-Berro et al. 
2004, Lee et al. 2004). We note, however, that these MACHO IMFs are more significantly 

enhanced with IM stars than the models considered here. While their goal was to generate a 
significant fraction of the Galactic Halo in white dwarfs, our intention is to produce enough 
Mg isotopes from IM stars to imitate a variation in a. Since this model quickly vanishes at 
early times due to the exponential decay factor in equation (8), it only can partially populate 
the Galactic Halo. An important consideration since the earlier work on the MACHO IMF 
models is that the estimated age of the universe has been reduced so that the white dwarfs 
would have had less time to cool. A smaller estimated age of the universe could be offset by 
a mass distribution that is peaked at a higher mass. 

Furthermore, MACHO IMF models are probably still difficult to reconcile with the C, N, 
and O constraint (Gibson & Mould 1997) and new observations of old white cool halo white 
dwarfs Flynn, Holopainen & Holmberg (2003). In spite of these difficulties with the halo 
MACHO models, we emphasize that the present models for Lyman-a systems are not affected 
by these constraints, both because the present models produce fewer white dwarfs than the 
MACHO models, and because the evolution of the Lyman-cc systems is not necessarily the 
same as that of the Galactic halo. 
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5.2. Observational Evidence of AGB-heavy IMF at low metallicity 

There is a large body of evidence suggesting a strong early contribution from IM stars. 
Indeed, ejecta from such AGB stars is needed to account for the observed abundances in 
both the low-metallicity stars in or near the Milky Way and those deduced in DLAs. Abia 
et al. (2001) give an overview of the abundances of the most extremely metal-poor stars 
in our galaxy and argue that an IMF peaked in the IM range of 3 — SMq is favored over 
an IMF populated by VMOs. This argument is based on the observed large number of 
low-metallicity halo stars with greatly enhanced C and N abundances. Intermediate-mass 
stars are particularly efficient at producing very large [C,N/Fe] ratios found in ~ 1/4 of the 
iron deficient stars (Norris, Beers, & Ryan 2000; Rossi, Beers, & Sneden 1999). The most 
likely source of these abundances is accretion from a binary AGB companion. Reinforcing 
this conclusion are the observations that s-process elements, especially Pb, appear even at 
the lowest of metaUicities (Aoki et al. 2001; Sivarani, Bonifacio, & Molaro 2003, Sivarani 
et al. 2004). Only in some cases can the s-process enrichment be attributed to Roche- lobe 
overflow from a companion. In others, an alternative explanation may be necessary. Because 
massive stars pollute the ISM before IM stars of the same age, the s-process cannot "outpace" 
the r-process except in several special circumstances: The r-process elements may not have 
mixed efficiently into the ISM and cluster; or the local first epoch of star formation may 
have been dominated by IM stars as described herein. 

It has also been argued that in order to account for the depletion of Deuterium, but still 
produce the observed white-dwarf population found the Galactic halo, a non-standard IMF 
is required (Fields et al. 1997; 2001). In particular, an early IMF dominated by IM stars has 
the advantage of matching these characteristics of galactic protodisks without the problem 
of metal over production from massive stars. In fact, our star formation rate, Eqs. (6) - (8), 
is based on the IMF proposed in Fields et al. (2001). 

There are two obvious constraints on the magnitude of the enhancement in the AGB 
IMF at low metallicity. These are carbon and nitrogen production, and the implied Type la 
supernova rate. An early IMF enhancement has an insignificant effect on the Type la rate 
because of the adopted minimum metaUicity requirement of Nomoto et al. (2003). By the 
time the ISM reaches the minimum metallicity, the IMF is dominated by the Salpeter IMF. 
An early IMF enhancement would, however, inevitably produce more carbon and nitrogen. 
Indeed, recent abundance determinations in low-metallicity systems have already shown a 
significant enhancement in carbon and nitrogen. 

Two such systems are the low-metallicity {Z ~ 0.004) globular clusters 47 Tucanae and 
M71, which show a strong pollution of AGB ejecta early in their evolution, yet maintain 
a constant iron abundance (Briley et al. 2003; Harbeck, Smith, & Grebel 2002). Many of 
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the observed stellar abundances show C depletion correlated with strong N enhancements. 
This trend is observed in both red giants and main-sequence turnoff stars, thereby making 
consistent internal contamination unlikely. Because globular clusters are sensitive to feedback 
from supernovae, we argue that the earliest star formation had to be comprised of mostly IM 
stars and not massive stars. Otherwise, the massive stars would contribute to the ISM first, 
and the C,N enhancements with respect to iron group would be diluted. Furthermore, SNe 
and the UV contamination by massive stars would suppress the star formation rate through 
feedback processes. We note that this model does not contradict the possibihty of a very 
early population of very massive stars which provide an effective prompt initial enrichment 
(Z ~ 10~^) and possibly reionize the Universe at redshifts z > 6. 

The Type la supernova rate does not depend upon our AGB enhanced IMF because of 
the adopted minimal metallicity constraint on the efficiency of binary accretion. Roughly 
2/3 of the iron produced comes from Type II and the rest is from Type la, as is consistently 
implemented in other chemical evolution models. Matteucci & Rccchi (2001) determined that 
a metallicity condition on Type la supernovae varies from the observations of the abundances 
in the Solar neighborhood. Allowing for Type la supernovae in our model at the earliest 
metallicities would shift the abundance evolution towards higher [Fe/H]. While the total 
abundance of Mg does not impact the QSO absorption spectra, a comparison of the degree 
to which an AGB enhanced IMF is distinguished from the constant IMF is detailed below. 
Both theoretical and observational studies of low-metallicity systems seem to point to an IMF 
dominated by AGB stars and their associated ejecta, at least in some systems. Competing 
hypotheses for the early IMF such as a constant IMF or VMOs have advantages in certain 
situations but often must also be fine-tuned to account for the full gamut of observations. 



6. Results 

We next probe the consequences of our adopted chemical evolution model on the pro- 
duction of the Mg isotopes and hence on the deduced shift Sa of the fine-structure constant. 
We compare these results to a "standard" model of chemical evolution (Timmes et al. 1995), 
which does not include the yields of AGB stars. This is listed as model in Table 1. We 
also compute the abundances of other elements and compare them to observations of DLA 
systems. It is particularly important to consider nitrogen and carbon abundances because 
these are also ejected from the IM stars responsible for Mg production. 

We will begin with the model used in Ashenfelter, Mathews & Olive (2004) (all models 
are summarized in Table 1). For this model (Model 1 in Table 1), we compute the abundances 
of N, N/Si, and C/0 as a function of [Fe/H], [Si/H], and [0/H] respectively. We then explore 
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the sensitivity of our results to our assumed star formation rate. Prom this study, we obtain 
a set of best-concordance parameters which reproduce the observed DLA abundances while 
still accounting for the data of Murphy et al. (2003a) (Model 2) or Chand et al. (2004) 
(Model 3) without varying the fine-structure constant. 
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Table 1: Model Parameters 

As shown in Ashenfelter, Mathews & Olive (2004), the birthrate function, Eqs. (6) - 
(8), provides for a large enhancement in the production of ^^Mg and ^^Mg due to the copious 
yields for these isotopes in IM stars at low metallicity. This result is shown in Pig. 4. The 
evolution of the heavy Mg isotopes can be explained as follows: Initially, the production 
of 25,26]y[g ^Yie ejecta of IM stars is delayed by their relatively long lifetime (compared 
to massive stars). Initial contributions to the chemical enrichment of the DLA interstellar 
medium comes from the most massive and shortest lived stars. In this model, the burst of 
IM stars begins to produce ^^'^^Mg at [Pe/H] > —2.5. At this stage, during the IM star- 
formation burst, ^^Mg and ^^Mg are efficiently produced relative to ^"^Mg as per the yields 
of Karakas & Lattanzio (2003). 

At higher metallicity, the ejecta from the standard population of (massive) stars (which 
is depleted in ^^'^^Mg) begins to dilute the ratio relative to ^"^Mg. This accounts for the decline 
from local maximum in ^^'^^Mg/^^Mg around [Pc/H] ^ — 1.5. At late times, the impact of 
the early generation of IM stars is diluted by the ejecta from subsequent generations of stars. 
The dashed curve excludes the burst of IM stars, while the dotted curve excludes the AGB 
yields as well as the IM component. This latter curve essentially reproduces the result of 
Timmes et al. (1995), which we refer to as Model 0. We note that the new AGB yields 
were also included in the chemical evolution model of Penner et al. (2003) who utihzed a 
normal stellar IMP. Their results are similar to that shown by the dashed curve in Pig. 4. 
While these results show higher abundances of ^^'^^Mg relative to ^^Mg than that given by 
the dotted curve, they are not high enough to account for the apparent variability in a. 

As discussed in section 4, the degree to which the heavy isotopes of Mg can effectively 
shift the value of the deduced fine-structure constant, depends sensitively on the ^^Mg/^^Mg 
ratio. This evolution of this ratio for Model 1 is shown in Pig. 5. As one can see, the 
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Fig. 4. — The chemical evolution of the ^^'^^Mg isotopes relative to ^^Mg. The solid curve is 
our result based on Eq. (6) using the AGB Mg yields of Karakas & Lattanzio (2003) and the 
parameters of Model 1. The dashed curve is the result obtained when the burst of IM stars 
is excluded while still including the AGB Mg yields. The dotted curve (Model 0) excludes 
both the AGB Mg yields and the enhanced IMF. The horizontal lines indicate the range for 
the ratio of ^^'^^Mg/^^Mg necessary to explain the shifts seen in the many-multiplet analysis 
as described in Section 4. 

26Mg/25Mg ratio remains slightly less than 1 for most of the metallicity range. This is due 
to the large enhancement of ^^Mg over ^^Mg in ~ 6Mq stars. The fact that this ratio is 
slightly less than Solar, requires a slightly larger excess of ^^'^^Mg relative to ^^Mg to explain 
the shift in a. 

The most vital comparison we must make involves the determination of the effective a 
variation in our model as a function of metallicity. Using the results of section 4, shown in 
Fig. 1, we can determine the shift in Sa based on the calculated Mg isotopic abundances. 
This result (for Model 1) is shown in Fig. 6. As one can see, the existence of the enhanced 
isotopic abundances of ^^'^^Mg at a metaUicity of [Fe/H] ~ — 1.5 implies no true variation 
in the fine-structure constant. DLA systems like the ones used to measure 6a/ a span a 
wide range of mctallicities, —1.75 < [Fe/H] < —0.75 with an average value of [Fc/H] ~ —1.1 
(Pettini 1999). By construction, this coincides with area of maximum effect from isotopic 
variations of Mg in our model. The horizontal lines show the shift in a needed to compensate 
for the variation deduced by Murphy et al. (2003a). 
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Fig. 5. — The chemical evolution of the ratio of ^^Mg/^^Mg. The curves correspond to the 
chemical evolution models as described in Fig. 4. The Solar ratio is 1.1. 



Since C and N are products of AGB stars, their observed abundances in DLA systems 
are capable of placing strong constraints on this particular chemical evolution model. In Fig. 
7, we show the N abundance produced in Model 1 with and without the AGB contribution. 
Comparing the results with observations of the nitrogen abundance in DLAs indicates that 
the AGB enhancement in Model 1 overproduces nitrogen. Fig 7 also shows that the standard 
model without an enhanced IMF (model 0) underproduces nitrogen. While this apparent 
underproduction may be the result of the stochastic nature of the the star formation process, 
in the context of our models at least some enhancement of the IMF for IM stars seems 
required by these data. The model parameters can in fact be adjusted to minimize the 
discrepancy with [N/H], while maintaining the requisite Mg isotopic enhancement to account 
for the variation in a. Indeed, the result shown in Fig. 6 is clearly dependent on our choice 
of the six basic model parameters listed in Table 1. In order to optimize the model with 
respect to both its ability to account for the apparent shift in a and to fit the abundances 
of the CNO elements, we vary each of the these parameters. The result of this optimization 
will be referred to as Model 2. The results of these variations are shown in Figs. 8 - 13, 
where we illustrate the influence of each parameter on the shift in 5a as well as their effect 
on the nitrogen abundance evolution as a function of [Fe/H]. In each of these figures, the 
solid curve corresponds to Model 1. 

In an independent high-resolution observation, the variation of the fine-structure con- 
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Fig. 6. — Shift in the fine- structure constant due to isotopic abundance variations in Mg. 
The shift in baja relative to an assumed Solar Mg composition is shown for both models 
with (solid) and without (dashed) an AGB-enhanced IMF. The shift needed to compensate 
the value ascertained by Murphy et al. (2003) is given by the horizontal dot-dashed line with 
the ±1(7 errors given by the horizontal dotted lines. 



stant was also measured by Chand ct al. (2004). They ascertained that baja = -0.06 ± 0.06 
assuming Solar Mg isotopic abundances. In order to completely account for this variation, 
a much smaller IMF enhancement would be needed. A model that fits their results also 
fits the nitrogen abundance well compared to a model with a standard IMF. We have also 
performed an optimized fit to their independent result. This is given as Model 3 in Table 1. 

To begin with, the star formation efficiency plays a major role in both the degree to 
which the heavy Mg isotopes arc produced as well as the degree to which nitrogen is produced 
in these models. The parameters for the efficiency of star formation and the infall time for the 
gas can be freely adjusted as well, since the galaxies corresponding to the DLAs encompass 
a wide range of morphologies and conditions. Timmes et al. (1995) determined that the 
star-formation efficiency parameter, esFi for the Milky Way in the local neighborhood was 
close to 3. Elliptical galaxies are generally regarded as being more efficient in forming 
stars than spirals, while smaller galaxies should be less efficient. Although the Mg isotopic 
ratio is largely insensitive to the star formation efficiency, it does affect the position of the 
peak with respect to [Fe/H]. A higher star-formation efficiency will produce more massive 
stars before the longer-lived IM stars recycle into the interstellar medium (ISM) as Fig. 8 
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Fig. 7. — Comparison of nitrogen evolution to observations in DLA systems. With an AGB 
enhancement (sohd hne), an excess of [N/H] is produced, while a model without an AGB 
enhancement (dotted line) underproduces nitrogen. These data were taken from various 
abundance measurement of individual DLA systems (Centurion et al 2003; Pettini et al 
2002; D'Odorico & Molaro 2004; Molaro et al. 2001; and Levshakov et al. 2002). 

illustrates. Relaxing the metallicity requirement for Type la supernovae would shift this plot 
and others to higher [Fe/H]. However, this adjustment can be compensated by reducing the 
star formation efficiency. 

In Fig. 8, we show the effect of varying esF on (^^Mg+^^Mg)/^'^Mg and [N/H] versus 
[Fe/H]. In addition, we see from Fig. 8b that indeed the trend in [N/H] abundances is better 
reproduced at low metallicity when esF is larger. 

As described earher, several model parameters affect the size of the IMF enhancement. 

The parameter A2 governs the weighting of the AGB enhancement with respect to the stan- 
dard Salpeter IMF component before the AGB enhancement has been exponentially reduced. 
As expected, the peak of the Mg isotopic ratio is linearly dependent on the parameter A2 
as seen in Fig. 9. In effect, the Mg isotope peak can be scaled with this parameter. The 
nitrogen overproduction is tempered when A2 is reduced. 

The time constants, ri and T2 also affect the total number of IM stars beyond the 
Salpeter IMF alone as shown in Figs. 10 and 11. The constant Ti determines the timescale 
for the onset of the normal component, and as one can see, has an important effect on 
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Fig. 8. — Dependence of the magnesium isotopic ratio (left) and the nitrogen abundances 
(right) on the coefficient of star formation efficiency, esF- All parameters were kept constant 
except esF: which had the values of 1 (dotted), 1.9 (solid), 3 (dot-dashed), 5 (dashed), 7 
(dot-dash-dashed), and 10 (dot-dot-dashed). 

the nitrogen abundance. At the same time, Ti causes the Mg isotope ratio to become less 
favorable toward accounting for the apparent variation in a. Fortunately, other parameters 
can boost the abundance of the heavy Mg isotopes without adversely affecting N/H. The time 
constant, T2 which determines how long the IM enhancement lasts is a prime example. By 
reducing Ti (0 corresponds to the standard component present at the onset of star formation) 
and increasing T2, we can maintain the high abundance of heavy Mg isotopes and at the same 
time reduce the nitrogen abundance. 

Other parameters do not change the total enhancement of IM stars; rather, they affect 
the distribution. The isotopic ratio is strongly dependent upon the center of the log-normal 
IMF enhancement, nic. Because the yield of the Mg isotopes is most efficient for progenitor 
stars with mass between 5— 6M0, the Mg peak will rise and fall depending upon the proximity 
of the peak in the IMF to this range. In Fig. 12a, the model parameters arc kept constant 
except nic- This figure shows that independent of the total size of the IM enhancement, the 
Mg isotopic ratio can be greatly enhanced by the choice of ttIc- The effect of the variation 
of TTic on the nitrogen abundance is shown in Fig. 12b. 

As Fig. 12a illustrates, the Mg ratio requires a value of rric around 5-6 Mq, whereas 
models with such high values for rric do more poorly for the evolution of N/H with respect to 
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Fig. 9. — Dependence of the magnesium isotopic ratio and nitrogen abundance on the 
coefficient of the IM component of the IMF, A2. All parameters were kept constant except 
A2, which had the values of 1 (dot-dot-dashed), 3 (dotted), 5 (dot-dashed), 8.9 (solid), 11 
(dashed), and 15 (dot-dash-dashed). 



the data. To further constrain the optimum value of m^, we can use a broader distribution 
of the AGB IMF component. However, the choice of ttIc does constrain the allowed range 
of the model parameter, cr, given in Adams & Laughhn (1996). In Fig. 13, we show the 
corresponding effects of the variation of the width a on the Mg isotopic ratio and nitrogen 
evolution. 

Finally, using our understanding of the available parameter space, the chemical evolution 
model was adjusted to obtain an optimal concordance with the measured abundances in 
DLA systems. While these systems exhibit a large degree of scatter in their respective 
abundances, a qualitative comparison can be performed. This comparison may also detail 
the consequences of adopting an IMF enhanced with IM stars, especially when compared to 
the standard Salpeter IMF model. Fig. 14 shows a comparison of calculated and observed 
[N/H] vs. [Fe/H] for the various chemical evolution models described herein. The optimized 
models, 2 & 3, provide a fair reproduction of the overall trend in the observed abundances. 

Fig. 15 shows the evolution of nitrogen relative to the alpha element silicon, which is 
dominantly produced in Type-II supernovae. This figure also illustrates the large abundance 
dispersions typical of DLA systems. The same models are compared to the [N/a] abundance 
as [a/H] evolves. The largest sample for Si abundances comes from the 20 measurements 
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Fig. 10. — Dependence of the magnesium isotopic ratio (left) and nitrogen abundance (right) 
on the turn-on time-scale of the standard component of the IMF, ti. All parameters were 
kept constant except ti, which had the values of (long-dashed), 0.05 (dot-dashed), 0.1 
(dotted), 0.2 (dot-dot-dashed), 0.5 (solid), and 0.7 (short-dashed). 



from Centurion et al. (2003). Although there is a great deal of scatter in these data, it 
has been interpreted that these abundances show a low dispersion [N/H]=-1.5 plateau and 
a second high dispersion plateau for [N/H]=-l. They contend that primary nitrogen from 
very massive objects would not be able to reproduce the lower plateau, while primary-yields 
from lower-mass stars do in fact produce this plateau. Both the AGB enhanced chemical 
evolution models and the standard IMF model show a plateau, however, the lower plateau 
is underestimated by 0.5 dex. This fact may indicate that the adopted WW95 primary 
yields underestimate nitrogen, or that an additional source of primary nitrogen has been 
overlooked. Another important consideration is that the Si yield depends sensitively on the 
adopted explosion energy and mass cutoff of Type 11 supernovae. The adopted choice of 
explosion energy in the no-metallicity model of WW95 leads to a dearth of Si produced. 
Consequently, this leads to a downward trend at the lowest values of [Si/H] in Fig. 15 
instead of the observed plateau. Centurion et al. (2003) point out that the yields of Meynet 
& Maeder (2002) match this plateau well. They further contend that the second plateau is 
due to the addition of secondary N from IM stars. While we do not observe a higher plateau 
in either the enhanced or standard IMF case, the enhanced IMF is more consistent with 
these abundances than the standard IMF. While both models suffer from a steep slope in 
between the two plateaus (the steep slope can explain the expeditious transition from the 
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Fig. 11. — Dependence of the magnesium isotopic ratio and nitrogen abundance on the decay 
time-scale of the IM component of the IMF, T2. All parameters were kept constant except 
T2, which had the values of 0.1 (dot-dot-dashed), 0.2 (solid), 0.3 (dot-dashed), 0.5 (dotted), 
and 0.7 (dashed). 



lower plateau to the higher one), the enhanced IMF model can be easily parametrized to 
account for the range of values of the upper plateau. The standard IMF is unable to account 
for the broad range of values in the upper plateau. However, one must keep in mind that 
the standard IMF does not have the parameter degrees of freedom that the enhanced IMF 
does. 

Since carbon is also a possible product of IM stars, we also illustrate the evolution of C 
relative to O (primarily from SNII) in Fig 16. Here we see that essentially all of the curves 
give a similar good fit to the observed abundance (Akerman et al. 2004) of metal-poor halo 
stars. Here we see that all of the models do equally well. The basic trend of these data 
is that at first the C/0 diminishes as supernovae produce more oxygen relative to carbon. 
Later for [Fe/H]> —0.5, C/0 increases due to the late-time ejection of carbon from low-mass 
AGB stars. The enhanced early IMF considered here has little effect on these curves, except 
for Model 1 in which substantial early carbon enrichment causes the C/0 to be a bit too 
high near [Fe/H] ~ —1.5. Nevertheless, our optimum models 2 & 3 do well in comparison 
to the data. We note that while it is reassuring that the models fit the data for local halo 
stars, chemical models for DLAs and the local neighborhood could be very different. 

Finally, in Figs. 17a and b, we summarize our estimated shift in deduced baja relative 
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Fig. 12. — Dependence of the Mg isotopic ratio (left) and nitrogen abundance (right) on 
the value of rric- All parameters were kept constant except rric, which had the values of 2.0 
(dot-dot-dashed), 4.0 (dotted), 5.0 (thick sohd), 5.5 (dashed), 5.75 (dot-dash-dashed), 6.0 
(dot-dashed), and 7.0 (thin sohd). 



to that obtained from a Solar Mg isotopic composition. It is clear that either of these models 
provide a possible alternative to the interpretation of a time varying fine-structure constant. 
Furthermore, the properties of these models of early AGB enhancements to the IMF are 
suggested by various abundance measurements in DLA systems and metal poor stars. In 
this interpretation, the MM method results have provided important new evidence into the 
star formation that occurred in the early epochs of DLAs. 

Before concluding, we make two final points. First, we emphasize that there is consid- 
erable scatter in the data. Scatter is found in both N/H data as well as the data from which 
Murphy et al. (2003a) and Chand et al. (2004) infer a value for Sa/a for each individual 
absorber. Our model can only be viewed as an average star formation history over many 
individual DLAs. Indeed, because of the stochastic nature of star formation, particularly at 
early times, we might expect large variations in the production of IM stars, and hence the 
production of the heavy Mg isotopes. For this reason, we believe that the MM method may 
provide a unique window to the star formation history in DLAs. Second, we also emphasize 
that, although the Chand et al. (2004) data are consistent with no variation in a, this is only 
so, because they chose Solar abundances for ^^'^^Mg. As one can see from Fig. 17b, Model 
3 which fits this data still requires a strong IM component in order to produce nearly Solar 



-3 -2.5 -2 -1.5 -1 -0.5 
[Fe/H] 



[Fe/H] 



Fig. 13. — Dependence of the magnesium isotopic ratio evolution on dimensionless width of 
the log-normal distribution, a. All other parameters were kept constant except cr, which had 
the values of 0.03 (dotted), 0.07 (solid), 0.1 (dot-dashed), 0.15 (dashed), and 0.2 (dot-dash 
-dashed) . 

isotopic abundances of Mg. 

7. Summary and Conclusions 

We have made a study of possible relations among stellar nucleosynthesis, the galactic 
chemical evolution of damped Lyman-a systems, and the apparent detections of a time- 
varying fine structure constant. In particular, we have explored the important effects of 
high temperature thermonuclear burning in low-metallicity AGB stars. Wc have shown that 
ejecta from these stars could have had a dominant effect on the early galactic chemical 
evolution of the crucial Mg isotopes in DLA systems. 

We have explored a variety of models in which the early initial mass function favors 
the formation of IM stars. Such an enhanced contribution from early IM stars allows for 
sufficient modification of the Mg isotope ratios to explain the many-multiplet results without 
a time- varying fine structure constant. Such a modified IMF may to some extent be a simple 
parametrization of uncertainties in theoretical estimates of ejected yields from AGB stars, 
but it is motivated by both theoretical and observational constraints. 
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Fig. 14. — Comparison of nitrogen evolution to observations in DLA systems. With Model 
1 (solid curve), a marginal excess is produced of [N/H]. The standard IMF corresponding to 
Model 1 (short dash curve) underproduces N these DLA systems. Models 2 (dot-dash) and 3 
(dotted) match the observed [N/H] considerably better than their standard IMF curve (long 
dash). The data was taken from individual abundance ratios of DLA systems (Centurion et 
al. 2003; Pettini et al. 2002; D'Odorico et al. 2002; Molaro et al. 2001; and Levshakov 2001). 

To compare with the MM method results, we have utilized an approximate treatment 
that qualitatively relates computed Mg isotopic abundances to the deduced fine structure 
constant for DLA systems in the redshift range of 0.5 < ^ < 1.8. Although this is only an 
approximation to the MM approach, we have shown that it reproduces the basic conclusions 
of detailed analysis (Murphy et al. 2003a). There is a real need to redo the MM method 
analysis in the context of evolving isotopic abundances as derived here. Incorporating isotopic 
variation would help to better quantify the need (or lack thereof) for IM stars and AGB 
nucleosynthesis in DLA systems. We hope that the present work will stimulate efforts along 
this line. 

In the context or our schematic analysis, we have explored a variety of chemical evolution 
models with an eye toward unraveling the time-varying alpha mystery while still satisfying 
the available constraints from observed elemental abundances in DLA systems. We have 
concentrated on the chemical evolution of N abundances, which are also produced in ACB 
stars. We also considered C as well as O and Si largely from Type II supernovae, and Fe from 
SNIa and SNII. We find that the observed high nitrogen abundances in DLA systems indeed 
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Fig. 15. — Plot of [N/Si] for [Si/H] for models with and without an AGB enhanced IMF. 
The data is compared to the DLA abundances of Centurion et al. (2003). Model 1 (solid), 
Model 2 (dot dash), and Model 3 (dotted) are compared to the standard (without an AGB 
enhanced) IMF models for 1 (short dash) and 2& 3 (long-dash). This data points illustrate 
two plateaus in the distribution of [N/Si] equal to around -1.5 and -0.9. 

confirm the need for enhanced ejecta from low-metallicity AGB stars. Even so, our previous 
model (Ashenfelter et al. 2004), which attempted to explain the MM results of (Murphy et 
al 2003a) tends to overproduce nitrogen and is therefore constrained by the observations. In 
this paper, however, we report on a parameter search which considers both data sets. We find 
a new optimum model (Model 2 in the present work) which simultaneously fits the observed 
N/H, C/0, and N/Si trends vs [Fe/H] while still eliminating the need for the time- varying 
fine structure constant as deduced from the Murphy et al. (2003a) data for systems in the 
redshift range 0.5 < z < 1.8. 

At the same time, we have also constructed a new model (Model 3) that can account for 
the results of the independent MM method analysis of Chand et al. (2004), which indicate 
smaller apparent variations in the fine structure constant. Even though these authors claim 
results that are consistent with no variation in a, this conclusion is based on the assumption 
that the Mg isotopic ratio is equal to the Solar one. As we have shown, in order to obtain a 
Solar isotope ratio at low metallicity, we must again rely on the role of IM stars and AGB 
nucleosynthesis . 

One conclusion of the present study is that important tests can be made of the hypoth- 
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Fig. 16. — Plot of [C/0] for [O/H] for models with and without an AGB enhanced IMF. 
The data is compared to the stellar abundances (Akerman et al.). Model 1 (solid), Model 2 
(long-dash), and Model 3 (dotted) are compared to the standard (without an AGB enhanced) 
IMF models for 1 (thin sohd) and 2& 3 (short-dash). 

esis that AGB nucleosynthesis can account for the apparent variation in the fine structure 
constant. The best measurement (though probably impossible) would be to directly detect 
Mg isotopic abundances from spectral lines. In our picture, the apparent variations in a 
should correlate directly with the fraction of heavy Mg isotopes. Conversely, if heavy mag- 
nesium abundances are significantly depleted relative to Solar, then the MM method results 
are actually understating the variation in the fine structure constant. Furthermore, if suffi- 
ciently precise data could be obtained to distinguish the ^^Mg and ^^Mg abundances, then 
large enhancements of ^^Mg observed in some systems could be indicative of Mg production 
that is specifically attributable to the Mg-Al cycle. If so, large enhancements of ^^Mg may 
also be anti-correlated with Al abundances. 

We further suggest that nitrogen (and/or carbon) abundances provide an easier test 
of the present hypothesis with regards to an enhanced IMF for low-metaUicity IM stars. 
Nitrogen should be measured and correlated with Sa/a in the same DLA systems to which 
the MM method is applied. A correlation of [N/H] with the largest variations in a would 
argue in favor of the present hypothesis. As another test of an enhanced IMF, the highest 
[Fe/H] or [Si/H] (Si and Fe are very correlated) in DLA systems should exhibit a significant 
over abundance of [N/a] if our IMF is correct. One caveat in using Si abundances is its 
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Fig. 17. — As in Fig. 6, the shift in the fine-structure constant due to isotopic abundance 
variations in Mg corresponding to Model 2 (left), the fit to the Murphy et al. (2003a) result 
and to Model 3 (right), the fit to Chand et al. (2004). 



dependence on the explosion conditions for the Type II supernovae. Future work regarding 
these tests should examine the consequences of adopting these different explosion criteria 
and compare with different yield models. 

Determining whether or not a metallicity condition exists for Type la supernovae is 
decisively related to an IMF enhanced with intermediate-mass stars. It may very well be 
that the reduced efficiency of Type la supernovae is ofi^set by the enhanced numbers of 
intermediate mass stars, which may satisfy the previous work of Matteucci & Recchi (2001) 
and the metallicity conditions of Nomoto et. al (2003). Comparison of abundances of N, 
C, and Mg with [0/H] may be able to test the viability of IM mass enhanced IMF models 
independent of Sn la rates. 

One should interpret the large scatter in the inferred variation in a, as well as the 
observed variation in element abundances like Mg, in the context of the stochastic nature of 
the star formation process at low metallicity. Our models simply represent a global average, 
and the contribution of IM stars may very well vary in individual DLAs. 

Clearly, more work needs to be done in constraining the effect of chemical evolution 
on the interpretations of DLA observations. Nevertheless, we have established that at least 
some fraction of the deviations in the fine-structure constant deduced from the MM method 
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could be due to chemical evolution effects. Among work that needs to be done, more stellar 
models of low-metallicity AGB stars over a broader mass and metallicity range arc needed to 
quantify the model predictions, particularly in the extrapolated mass range of 7 < M < 12 
Mq between the AGB and Type II regimes. In particular, it would be very useful to have 
have a full complement of stellar yields (including C, N and Mg) in the IM range derived 
from a self-consistent set of stellar models. Out of necessity, we supplemented the Mg yields 
of Karakas & Lattanzio (2003) with the Padova CNO yield models; however, a self-consistent 
model would more accurately quantify the correlation of Mg from AGB sources to N. 

At the same time, more observations of abundances in DLA systems arc needed. Most 
importantly, the DLAs that are used in the MM-mcthod should have their associated abun- 
dances quantified in order to see if there are correlations between the apparent variations 
in a and various abundances. Additionally, more atomic physics work needs to be done to 
more accurately quantify the possible isotopic shifts in the absorption lines (particularly for 
the high-z data). 

The MM method has presented a very important question as to whether the fine struc- 
ture constant varies with time. It is hoped that the present study will stimulate further 
efforts along all of the above lines with a goal of clarifying this important physical question. 

If the fine structure constant docs vary in time or space, it provides an important window 
into the physics beyond the standard model. As such, the chemical evolution efi^ccts de- 
scribed herein should be carefully quantified to reduce the systematic uncertainties in the 
deduced result. Even if our chemical-evolution interpretation of the MM results prove to 
be verified, the MM method will have provided valuable insight into the mysteries of early 
cosmic star formation and galaxy evolution. 
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